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The Botvina and Mishustin's comment addresses a de- 
tail (minor subject to the aim and topic of the paper, 
i.e. the Coulomb effects on the liquid-gas phase transi- 
tion) of the recent paper namely a free volume ap- 
proach aimed to allow some qualitative estimation of the 
system's critical point. In their view this distorts the be- 
haviour of the pressure (P) versus volume (V) curves due 
to the "non-conservation of the center of mass (cm.) of 
the system" implied by the employed parametrization. 

While we agree with the evident idea that, apart from 
the case of the large systems, the center of mass of a small 
system should be conserved in thermodynamical calcula- 
tions in order to obtain accurate results, one should re- 
mark that the Botvina and Mishustin's argumentation is 
based on a specific view on freeze-out: non-overlapping 
spherical fragments placed in a spherical "recipient" of 
volume V . The shape of the fragments is therefore fixed 
for each fragmentation event. This is, indeed, a sce- 
nario widely used in statistical models from the beginning 
of the nuclear multifragmentation studies. Within such 
a scenario, the Botvina and Mishustin's comments are 
surely correct and they apply to most of the statistical 
calculations performed so far (MMMC model Q , Koonin 
and Randrup model Das Gupta and Mekjian calcu- 
lations 0, and even recent calculations performed with 
the SMM model etc.). Following this hypothesis, at 
a given temperature, P will converge to a finite value 
when V will approach the volume of the source nucleus 
at normal density (Vq) (see e.g. 6]). 

However, dynamical calculations clearly show that 
fragments are not at all spherical at break-up but have 
various shapes variating from event to event (see e.g. 
It means that the spherical fragment hypothesis is actu- 
ally away from reality and this has to be reflected in all 
thermodynamical observables. In this case, the system's 
position coordinates (entering the integration which leads 
to the free volume term in the partition function) will be 
shared between the positions of each fragment's cm. and 
the additional coordinates corresponding to each frag- 
ment deformation. (In other words there are more than 

— 1 contributing degrees of freedom.) Obviously, this 
detail will raise the power of Vf to a value larger than the 
A^ — 1 responsible for the finite value of P when V ^Vq, 
and therefore P will diverge near this limit leading to 
the well known van der Waals behavior (see for example 



the nuclear phase diagrams from Ref. [8j). This very 
simple reasoning shows that the conclusions of Botvina 
and Mishustin are actually in disagreement with what 
one would expect from a real nuclear system. 

But let us now discuss the free volume approach we 
used Ref. This approach was employed in Ref. 
with the intention of reaching densities p po unaf- 
fordable with a spherical fragment scenario. Within this 
scenario it is assumed that each fragment is blocking for 
itself and the rest of the fragments its own volume. Ob- 
viously, at p — > po such a formula can only be justified 
by a deformed fragment scenario. And this differs a lot 
from what Botvina and Mishustin claim to understand 
from our paper: the removal of the physical assumption 
that fragments don't overlap. The latter assumption is 
indeed unphysical while the deformed fragment scenario 
(actually used in our calculations) is not. 

Of course, the deformed fragment scenario deserves a 
more refined treatment than our calculations from Ref. 

a proper account for the fragments' deformation de- 
gree of freedom, a sharp conservation of the center of 
mass, an accurate treatment of the binding energies, etc. 
Definitely, these aspects have to be addressed in the sta- 
tistical models for a better description of the real multi- 
fragmentation. 

Summarizing, the Botvina and Mishustin's P{V)t 
curves correspond to a simplified scenario for the freeze- 
out: nonoverlapping spherical fragments inside a spher- 
ical freeze-out recipient, which, as explained before dif- 
fers from the scenario corresponding to our free volume 
parametrization: non-overlapping fragments with vari- 
ous shapes at freeze-out. The latter scenario is the re- 
alistic one, being in agreement with the results of the 
various dynamical simulations of nuclear multifragmen- 
tation and, moreover, at p ^ po the fragments have to be 
deformed in order to fit the recipient. As explained be- 
fore, within this scenario the pressure has to diverge while 
V approaches Vq (fact in agreement with our results and 
in disagreement with Botvina and Mishustin ones) since 
the power of Vf will necessarily be greater than A^ — 1. 
More refined treatments than ours are surely necessary 
for obtaining very accurate results, but the qualitative 
character of our evaluation was already stated in Ref. 
Q. It results that, in spite of the roughness (commented 
above) of our approximation, the Botvina and Mishus- 
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tin's result [on the P{V)t curves] is in fact the "mislead- 
ing" one since it differs essentially from what is expected 
from real systems. 
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